Jammed disks in narrow channel: criticality and ordering tendencies 
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A system of identical disks is confined to a narrow channel, closed off at one end by a stopper and 
at the other end by a piston. All surfaces are hard and frictionless. A uniform gravitational field is 
directed parallel to the plane of the disks and perpendicular to the axis of the channel. We employ 
a method of configurational statistics that interprets jammed states as configurations of floating 
particles with structure. The particles interlink according to set rules. The two jammed microstates 
with smallest volume act as pseudo- vacuum. The placement of particles is subject to a generalized 
Pauli principle. Jammed macrostates are generated by random agitations and specified by two 
control variables. One is a measure of the intensity of random agitations at given pressure. The 
other is a measure of the change in gravitational potential energy in units of compression work when 
one particle is excited. In this two-dimensional space of variables there exists a critical point. The 
jammed macrostate realized at the critical point depends on the path of approach. We describe all 
jammed macrostates by volume and entropy. Both are functions of the average population densities 
of particles. Approaching the critical point in an extended space of control variables generates two 
types of jammed macrostates: states with random heterogeneities in mass density and states with 
domains of uniform mass density. 



I. INTRODUCTION 

Granular matter consists of particles with sizes too big 
for thermal fluctuations to have a significant impact. The 
spatial configurations of grains depends on many factors 
including their sizes and shapes, interactions between 
them, interactions with external fields, and a protocol 
for producing them in a controllable way [TpE]. Models 
of granular matter tend to omit non-essential attributes 
of grains in efforts to gain more profound knowledge of 
bulk behavior. 

In this spirit, jammed states of rigid objects with no 
interactions other than hard-core repulsion have become 
the focus of numerous studies. It has become common to 
describe jammed states of granular matter in the lan- 
guage of equilibrium statistical mechanics albeit with 
provisos. We shall use the term configurational statis- 
tics for this framework of analysis [1HS]. Jammed states 
are frozen, hence time averages useless. Nevertheless, 
jammed macrostates are postulated to exist, to be sys- 
tematically reproducible by means of some protocol of 
random agitations, and to be describable by entropies 
and averaged densities in spatial regions of macroscopic 
size. 

To the many interesting themes surrounding jammed 
granular matter belong the existence and nature of crit- 
ical singularities and the onset of ordering [7]. In this 
work we use a very simple scenario that enables us to 
investigate this theme with some degree of rigor. 

Inspired by prior work [SMlO] , we consider a long, nar- 
row channel of width H with the axis in a horizontal 
direction, a stopper at one end, and a piston at the other 
end as shown in Fig. [T] The channel contains N disks 
of mass m and diameter a. All surfaces are rigid and 
frictionless. Jamming requires three points of contact on 
each disk that are not on the same semicircle. A uniform 
gravitational field g in the direction shown is present. 
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FIG. 1: Jammed microstate of disks of mass m and diameter a 
in a channel of width H with rigid walls, stopper, and piston. 
For 1 < H/a < 1 + \/3/4 all configurations are sequences of 
four distinct interlinking two-disk tiles. In our methodology, 
two tiles (marked v) are elements of pseudo- vacuum and the 
other two (marked 1,2) are particles. 



In channels with 1 < H/a < 1 + y3/4, every jammed 
disk (except the two outermost) touches one wall and 
two adjacent disks. Every possible jammed microstate is 
a sequence of interlinking two-disk tiles. The four dis- 
tinct tiles are identified in Fig. [TJ The two most compact 
states consists of tiles v in alternating sequence [TT] . The 
least compact states are sequences [vlv2vlv2 • ■ • v] and 
[v2vlv2vl • ■ ■ v]. Only tiles v can follow each other di- 
rectly. Any tile 1 (2) is separated from the next tile 1 
(2) by at least two tiles v. A tile 1 is separated from the 
nearest tile 2 by at least one tile v. 

In channels with 1 + -\/3/4 < H/a < 2, there exist 
32 tiles interlinking according to more complicated rules 
[3]. If H/a > 2 disks can pass one another with the con- 
sequence that the jamming condition becomes nonlocal. 
Here we assume that H/a < 1 + a/3/4. 

Jammed microstates in this system are countable. 
Each has a well defined volume. Many have the same 
volume. They represent one jammed macrostate in the 
sense of an ensemble average. The multiplicity of mi- 
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crostates with given volume determines the configura- 
tional entropy of that macrostate. In the limit N — > oo 
the entropy per disk, S = S/N becomes a smooth func- 
tion of reduced excess volume, V = (V — Vq)/N, where 
Vq is the volume of the most compact jammed state. 

The inverse slope of the configurational entropy curve 
is an intensive thermodynamic variable named com- 
pactivity [2]: 



X 



dS 
dV 



(1) 



This definition works well when a jammed macrostate is 
specified by a single intensive variable. Typical for such 
cases is that S increases monotonically between V = 
and some value V — V^. The function S(V) is concave 
along that stretch with infinite initial slope and zero fi- 
nal slope. The compactivity thus varies between X = 
at V — and X = oo at V — Voo- It is a measure 
for the intensity by which the system is randomly agi- 
tated at given pressure of the piston to produce a specific 
jammed macrostate. A dosage of intense random agita- 
tions produces a highly compactifiable state with much 
excess volume. Jolting the system with less intensity pro- 
duces a more compact and less compactifiable jammed 
macrostate. 

Between the most disordered jammed macrostate at 

V = Voo and the least compact jammed macrostate at 

V = V rn ax the function S(V) stays concave and thus de- 
creases. These jammed macrostates have negative com- 
pactivity. They do exist but cannot, in general, be gen- 
erated by random agitations. 

This scenario is borne out if we set g = in our sys- 
tem or tilt the channel into a horizontal plane. The 
configurational entropy of this case was determined in 
[5]. The gravitational field g complicates matters signif- 
icantly. Jammed macrostates are now specified by two 
independent intensive variables. Microstates with equal 
volume do no longer have equal statistical weight on ac- 
count of their different gravitational potential ene rgy . 

The methodology is further developed (Sec. Eb. S 
turns out to be a multiple- valued function of V (Sec~ [lll[ ). 
In the space of the two independent intensive variables 



there exists a critical point (Sec. IV I. Ordering tendencies 



at criticality are explored in a space of three independent 
intensive variables (Sec.|v]). More complex scenarios now 
appear within reach of rigorous analysis (Sec. VI). 



II. METHODOLOGY 

Here we adapt a method of statistical mechanical anal- 
ysis based on fractional exclusion statistics to jammed 
states of granular matter. This approach was originally 
developed for quantum many-body systems [T^HT5] and 
was recently extended to classical particles with shapes 
[T6HT9] . Whereas microscopic degrees of freedom are sub- 
ject to thermal fluctuations, the disks in the narrow chan- 
nel need to be subjected to random agitations to produce 



a comparable effect. In the context of this study the 
disks are not the particles themselves. The particles (or 
quasiparticles) considered here are the tiles 1 and 2. The 
method to be described has a combinatorial part and 
a statistical mechanical part. We begin with a discus- 
sion of the variables that we employ to describe jammed 
macrostates. 



A. Energetics 

The volume of each element of pseudo- vacuum as in- 
dicated by a bracket in Fig. [l]is a-y/l - (H/a-l) 2 and 
the volume of a particle from either species is a. Hence 
the activation of each particle increases the volume by a 
fraction of the disk diameter, 



vi = v 2 



qa, 



1 - - (H/a If 



(2) 



The activation of a particle 1 or 2 adds energy. One disk 
moves up or down and the piston moves out: 



ei,2= W±7; r y = mg(H-<r). 



(3) 



We shall see that the configurational entropy, S(V), 
now depends on a parameter that contains the ratio p/g- 
Jammed macrostates are specified by two independent 
intensive variables. We get to them naturally, in the 
framework of configurational statistics, if we start from 
the general thermodynamic relation (dS /dV)jj =p/T. 

For g = this is the inverse of compactivity |l]), the 
only intensive variable needed. In a thermal system, the 
dimensionless quantity, 



p ~ k B r 



(4) 



compares a unit of thermal energy with a unit of com- 
pression work. For the athermal system of jammed disks 
this translates into a measure of the intensity of random 
agitations performed at a given pressure of the piston. 
The jammed macrostate resulting from random agita- 
tions only depends on that ratio. 

For j / a second dimensionless quantity comes into 
play. In a thermal system the quantity, 



K =^L_ 

9 ~ fc R T' 



(5) 



compares the same unit of thermal energy with a con- 
venient unit of gravitational potential energy. For the 
jammed disks this translates into a different measure of 
intensity of random agitations, performed in a given grav- 
itational field. Configurational statistics postulates that 
jammed macrostates are reproducibly characterized by 
K p and K g . 

The thermodynamic variables T and p per se have no 
meaning for the description of jammed macrostates. Nev- 
ertheless, in an experiment, random agitations are per- 
formed at some piston pressure. The resulting jammed 
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macrostate depends on the intensity of random agitations 
and on how that pressure compares with the weight of the 
disks. The entropy S and the excess volume V are func- 
tions of K p and K g , but the shape of the function S(V) 
inferred from them only depends on the experimentally 
controllable parameter, 



K9 



7 
pqa' 



(6) 



B. Combinatorics 

For the combinatorial analysis we use the template de- 
signed for statistically interacting particles familiar from 
previous work (adapted to open boundary conditions) 

[nun]: 



M 



d m + N m - 1 



W({N m })=n pv ]l l" m >£™ " J , (7a) 

m— 1 ^ ' 

M 

dm A m ^ ^ 9mm' (Nm' ^mm')' (7b) 



where the (real, positive) w m are the solutions of the 
coupled nonlinear algebraic equations, 



A I 



= (l+W m ) II C 1 + «£')"'"'"'> (10) 



and where K m = ^m/ksT. The average population den- 
sities, N m = (N m )/N, of particles from each species are 
derived from the coupled linear equations, 



M 



IV, 



iNm ^ Qmm' Nm' — A m s 



(11) 



where A m = A m /N. For jammed-disk macrostates the 
K m are linear combinations of Q and (J5j|. The con- 
figurational entropy, S(V), follows directly from a scaled 
version of the entropy expression ([8| in combination with 
the volume expression [2D] . 



M 



V({Nm}) = 



(12) 



III. CONFIGURATIONAL ENTROPY 



This multiplicity expression yields the number of mi- 
crostates with given particle content. The degeneracy 
of the pseudo- vacuum is encoded in n pv . The A m are ca- 
pacity constants and the <? mm < are statistical interaction 
coefficients. 

All particles of a given species considered here will have 
the same volume and energy no matter where they are 
placed. Therefore, all microstates with given particle 
content constitute a macrostate in the sense discussed 
earlier. The entropy of a macrostate as derived from the 
multiplicity expression ^ for N m >• 1 via S — k B \nW 
reads [HUH]: 



M 

S{{N m }) = k B [(Nm + Y m ) In (jV m + Y m ) 

(8a) 

M 

(8b) 



Y m Am ^ Qmm'N^ 



m' — l 



All jammed microstates such as the one shown in Fig. [I] 
contain particles from M = 2 species. The specifications 
needed for the statistical mechanical analysis are listed 
in Table [I] In the taxonomy of [17] both species are com- 
pacts. 



TABLE I: Specifications of the particles represented by tiles 
1 and 2: species, volume and energy (relative to vacuum), 
capacity constant (left), and interaction coefficients (right). 
The elements of vacuum have (absolute) volume (1 — q)a and 
energy p(l - q)a. 
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A. Zero gravity 



C. Statistical mechanics 

The statistical mechanical analysis for a thermal sys- 
tem |T3l [14] produces, for the grand partition function, 
the result [T7] 



A I 



(9) 



In the special case g = 0, particles 1 and 2 have equal 
energy. No distinction is necessary. They can be merged 
into a single species I with specifications A\ = N — 3 
and gn = 2. It is a type-1 merger in the classification 
of |19j . The configurational entropy follows directly from 
^ with these specifications: 

S = k B [(l-Ni) ln(l -Nf)-Ni In Aj 

-(1-2^)1^1-2^)], (13) 
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where V = gerAj. This result was first derived in [5] by 
a different method. In this case we do not have to solve 



Equations ( 10 1 in simplified form read 



(10 1 and (11) to get to S(V). Those solutions are 



1 + \J\ + 4e-- Rr i 



Ni = 



(14) 



where K\ = K p from Q. The entropy inferred from Q 
reads 



5' 



In (1 + 10= 



Ki 



(15) 



Expressions ([14| and ( |15| are consistent with (13). All 
the attributes of S(V) described in Sec.[l]are here realized 

m 



B. Nonzero gravity 

For g^O the entropy depends on two independent 
variables. The function S{N%,N2) inferred from (|8| with 
the A m and g mm ' from Table [I] produces an entropy land- 
scape as shown in Fig. [2j It is of quadrilateral shape with 
a smooth maximum in the center and zeros at the cor- 
ners. All jammed states on a line N% + N2 — const have 
the same volume bu t di fferent gravitational potential en- 
ergies. Expression |l3| ) is recovered by S(Ni/2,N%/2). 
The system has the highest capacity for particles if both 
species are present in equal numbers. 
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FIG. 2: Entropy per disk, S/ks, versus population densities 
Ni,Na. The contours are at 0.049£, t = 1, . . . , 9. The seven 
curves represent particle population densities of macrostates 
with r = 0, 0.5, 0.75, 1, 1.25, 1.5, 3 (left to right). 
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W\W2 



(1 + wi){\ + w 2 ) 



Wi 

U>2 



(16) 



They reduce to a 4 th -order polynomial equation with one 
physically relevant, real solution w m (K p , K g ), m = 1,2. 
The particle population densities inferred from (11) arc 



w 2 



D 



K n 2 = Wi 



1 



D = 2wiW2 + 3(u>i 



D 1 
w 2 ) +4. 



(17a) 
(17b) 



Jammed macrostates generated by random agitations 
fall onto lines parametrized by T. Several such lines are 
shown in Fig.[2j All lines begin at the jammed macrostate 
with maximum entropy and equal population densities. 
With K p increasing from zero the lines fan out. On the 
line r = the two population densities are being depleted 
at the same rate. Higher values of T have the effect that 
particles 1 are depleted more rapidly. The population of 
particles 2 may decrease or increase as the curves show. 

At r = 1 particles 2 have zero activation energy. Their 
density increases moderately. This is an entropic effect. 
As the population of particles 1 (with positive activation 
energy) decreases, the entropy is maximized by an in- 
creasing population of particles 2. The (nearly straight) 
line pertaining to this case terminates at the point of 
maximum entropy for states with N\ = 0. The configu- 
rational entropy curves, S(V), as inferred from (|16j) and 
(17) in combination with ((8|) are shown in Fig. |3| for the 



same values of T. 

The concave curve on the upper left represents the 



zero-gravity case discussed in Sec. Ill A and previously 
solved in [8 . The inverse slope of that curve is the 
compactivity ([I]), a quantity that varies between zero in 
the most compact jammed macrostate and infinity in the 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
V/qo- 



FIG. 3: Entropy per disk versus excess volume at T = 
0,0.5,0.75,1,1.25,1.5,3 (left to right). The dashed arch de- 
limits a region of critical macrostates. 
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macrostate with the highest possible entropy. All curves 
for nonzero gravity start from that highest-entropy point 
and fan out in different directions. 

At r < 1 both species of particles have positive activa- 
tion energies. The physics does not change qualitatively 
from the zero-gravity case. The gravitational field sup- 
presses the population density of particles 1 and enhances 
that of particles 2. The net effect is a modest reduction 
of entropy (of mixing) at given volume. All curves end 
in the most compact jammed state. 

At r > 1, the particles from species 2 have negative 
activation energies. The jammed state produced by ran- 
dom agitations in the low-intensity limit is now quali- 
tatively different. It is a periodic array of close-packed 
particles 2. It has volume V/qa = |, which is larger 
than the (r-independent) jammed state produced in the 
high-intensity limit of random agitations. 

Increasing the intensity of random agitation always 
increases the entropy. However, at T > 1, it initially 
decreases the volume of the states thus produced. The 
decreasing V is associated with compression work used 
for weight lifting. The annihilation of particles 2 creates 
disorder and thus increases S. At higher intensities the 
jammed states have larger volume again. Here particles 
1, which have high energy, are created in significant num- 
bers. 

For the borderline case T = 1 between the two regimes, 
the particles from species 2 have zero activation energy. 
Gravitational energy and expansion energy are in bal- 
ance. Random agitations do not produce ordering in the 
low-intensity limit. The shapes of the curves in Fig. [3] 
signal the presence of a critical singularity. 



Further critical macrostates are represented by points 
on the (dashed) arching line in Fig. [3] and by points in the 
space underneath. The latter will be discussed in Sec.|V] 
The former are realized in a family of pathways toward 
the critical singularity specified by a single parameter, 



A = K P -K g . 
The critical macrostates are determined from 



(22) 



Wl 



00, 



,2A 



I + W2' 



-00 < A < +00 (23) 



and the particle population densities (111 become 
Aq =0, N 2 = 



1 



3 + 2w 2 



(24) 



The point A = recovers the case T = 1 discussed pre- 
viously. 

Only particles 2 are present. They are distributed ran- 
domly along the channel. The volume and the entropy 
depend on the population density N 2 : 



V/qa = N 2 

S l k B = ^ 

1 - 



hi 

3A 2 
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(25a) 



(25b) 



Complete ordering exists if there are no particles or if the 
particles are close packed. The function S(V) has zero 
slope at the point (21 1 and infinite slope at S = 0: 



IV. CRITIC ALITY 

In the framework of configurational statistics it is le- 
gitimate to use the term criticality more loosely than is 
common in the theory of (thermal) phase transitions. In 
our system the critical singularity is associated with the 
combined limit 



K„ 



k„ -> 00, r ->• 1. 



(18) 



The jammed macrostate realized at the critical singu- 
larity depends on how the singularity is approached. In 
Scc. |III| we have considered one particular approach, with 
r = 1, implying that w± 00 and 



w 2 



(9 ~ V09) 1/3 + (9 + VG9) 1/3 



1.3247..., (19) 



2V332/3 

the physically relevant solution of the cubic equation, 



w 2 



1 = 0. 



(20) 



The associated values for excess volume and entropy are 

S/k B ~* In w 2 =0.28119..., (21a) 
V/qa ~* (3 + 2w 2 ) _1 = 0.17700 .... (21b) 



dS 
dV 



- In V, 
-|ln( 



V ->• 
V), v^\ 



(26) 



All critical macrostates are averages of the same sub- 
set of jammed microstates, but weighted differently. The 
macrostate at the top of the arch has all microstates 
from that subset weighted equally. The macrostates else- 
where on the arch result from averages with one bias in 
the weighting. The parameter A controls the population 
density N 2 and thus the volume. All microstates with a 
given number of particles 2 remain weighted equally. The 
macrostates inside the arch result from averages with two 
biases in the weighting. 



V. ORDERING TENDENCIES 

We choose a second bias in the form of a weak inter- 
action force that either enhances or suppresses hetero- 
geneities in mass density along the channel. It controls 
the degree of clustering of particles 2. New pathways to 
the critical point generate macrostates with domains of 
uniform mass density. No symmetry-breaking field needs 
to be introduced to make that happen. 
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We extend the set of particle species from two to three. 
We split one species into two: every compact 2 in the 
system becomes either a host 2' or a tag 2". The particle 
2 closest to the piston is a host 2'. Any other particle 2 is 
a tag 2" if it follows the preceding particle 2 in one of the 
two shortest tile sequences: 2vv2 or 2vlv2. Otherwise it 
is a host again. 

In Fig. [I] we have one host 2' followed by a tag 2" 
(left to right). We do not split compacts 1 because they 
become depleted as criticality is approached. The spec- 
ifications of the three species, compact 1, host 2', and 
tag 2", are compiled in Table [Tl| Anghel's rules [TS] are 
satisfied. 



TABLE II: Specifications of three species of particles: species, 
volume and energy (relative to vacuum), capacity constant 
(left), and statistical interaction coefficients (right). 
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Raising (lowering) the activation energy of the tag rel- 
ative to that of the host amounts to a repulsive (attrac- 
tive) short-range force between particles 2. Such a force 
suppresses (enhances) heterogeneities in mass density. It 
opens up new pathways to the critical singularity. In 
the critical macrostates thus generated the microstates 
with a given number of particles 2 are no longer weighted 
equally. The weighting now discriminates between hosts 
and tags. 

The evidence for the formation of domains presented 
here is indirect, as found in the entropy |22j . Assembling 
the function S(V) in the extended parameter space starts 
from the extended Eqs. (16), 



„2K« 



w\wl,(l + W 2 ") 

(1 + Wl)(l +W 2 ,) 2 W 2 ii ' 

Wl(l +w 2 >)w 2 " $ _ 
1«|,(1 + W 2 ") 



w 2 > 

W 2 n 



(27) 



where $ = 0/fcsT is an amendment to Q and The 
population densities of compacts 1, hosts 2', and tags 2" 
are 

1 + 2w 2 n + W 2 >W 2 n 
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(28) 



7w 2 " + 5wiw 2 » 
- 3w 2 'iv 2 " + 2wiw 2 >w 2 > 



We now have a two-parameter family of pathways to- 
ward the critical singularity with critical macrostates de- 



Wi = 0, 



„2A 



W%{1 + w 2 ») 2 



w 2 < 



(1 + W 2 >) 3 W 2 » ' W 2 n 

oo < A < +oo, -co < $ < +oo, 



(29) 



and with critical population densities, 



Nl = 0, N 2 >= W 2 nN 2 „ 



W 2 i, 



3 + 5u>2" + 2w 2 >w 2 i. 



(30) 



The function S(N 2 i , N 2 ") inferred from ^ is shown in 
Fig. [4j Every point in this entropy landscape represents 
a macrostate at criticality. All critical macrostates have 
pqa = 7 and <f> = 0. 
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FIG. 4: Entropy per disk, S/ks, versus population densities 
N2",N 2 ' of particles from the two species 2" and 2'. The 
contours are at 0.02 81, £ = 1,...,9. The solid and dashed 
curves represent (301 for <f> = 0, ±1.5, ±3.0, ±5.0 in ascending 



order from top to bottom. 

In macrostates pertaining to $ = 0, particles 2 form 
random clusters. They are the states of highest entropy 
at given N 2 > ± N 2 n , represented by the dashed curve in 
Fig. [3] This particular mix of hosts and tags, 



Nv = 



1 



N 2 „(8 + N 2 „)-5N 2 



(31) 



produces a characteristic texture of heterogeneity in mass 
density. Critical macrostates on either side of the dashed 
curve have mass heterogeneities with textures that reflect 
a higher degree of ordering. 

In macrostates generated on pathways with $ > 
tags are more abundant and hosts less abundant than 
in a mixture of randomly placed particles 2. The bottom 



three solid curves are critical population densities ( 30 1 at 
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fixed $ > 0. Clusters of hosts with tags are domains of 
uniform, low mass density, phase separated from domains 
of uniform, high mass density. The latter are clusters of 
tiles v. 

These domains suppress the entropy significantly rel- 
ative to the entropy of random clusters of particles 2. 
The entropy reduction is demonstrated by the contours 
in Fig [4] and, more compellingly, in the configurational 
entropy curves, S(V), shown in Fig. ^a). The effect is 
most significant at intermediate volume where large clus- 
ters of tiles 2 or tiles v are least likely to be realized by 
chance. 

A different ordering tendency is manifest in critical 
macrostates generated along pathways with $ < 0. More 
hosts and fewer tags are present than in a randomly 
placed mixture. Critical population densities (30) at 



fixed $ < are represented by the top three solid curves 
in Fig. [3] These macrostates exhibit various degrees of 
dispersal of particles 2. The dispersal flattens out hetero- 
geneities in mass density and thus reduces the entropy. 

This ordering tendency manifests itself differently in 
the configurational entropy as shown in Fig. Rub). Dis- 
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persal of particles 2 only produces ordering tendencies 
in conjunction with spatial constraints. The entropy re- 
duction remains insignificant until crowding of particles 
2 becomes an issue. The capacity of the channel for hosts 
2' is only 60% of its capacity for particles 2. Hence, hosts 
alone reach saturation at V /qa — | in an ordered state. 
At larger V, the entropy rises again on account of tags 
attached to some hosts. 
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VI. DISCUSSION AND OUTLOOK 

The methodology of this work builds on the founda- 
tions of configurational statistics. Chief among the as- 
sumptions is the existence of jammed macrostates that 
are systematically reproducible by some protocol of ran- 
dom agitations of given intensity and that can be de- 
scribed by few control variables. 

Jammed macrostates of the system under investigation 
here are parametrized by two intensive variables reflect- 
ing the intensity of the random agitations in two distinct 
energy units of relevance. One unit represents the com- 
pression work of the piston when it is moved a fraction of 
one disk diameter. The other unit represents the weight 
of one disk lifted across the channel. 

A protocol of random agitations that fits the bill of con- 
figurational statistics for the jammed disks in the narrow 
channel is not guaranteed to exist. If it does exist it may 
be difficult to implement. It is no secret that this is the 
Achilles' heel of configurational statistics. 

The detailed predictions of this work, which follow rig- 
orously from the assumptions underlying configurational 
statistics, present a situation where these very assump- 
tios can be put to the test of simulations and experiments. 
The simplicity of the model in conjunction with the com- 
plexity of its behavior is well suited for that purpose. Any 
instances of verification or falsification will shed light on 
the important question of how accurately a protocol of 
random agitations applied to jammed granular matter 
can mimic the effects of thermal fluctuations. 

The methodology which combines combinatorial anal- 
ysis and statistical mechanical analysis in an adapta- 
tion to jammed granular matter as presented here is 
amenable to further development in several directions: 
(i) the regime 1 + -\/3/4 < H/ a < 2 of 32 tiles as identi- 
fied in [5], (ii) mixtures of disks with different radii and 
masses including rattlers [23], (iii) spatial correlations of 
heterogeneities in mass density, and (iv) channels with 
the axis oriented vertically. 



FIG. 5: Entropy per disk versus excess volume at critical- 
ity. The five curves, inferred from (291, represent, from top 



to bottom, critical macrostates characterized (a) by $ 
0, 1.5, 3, 5, 10. and (b) by $ = 0, -1.5, -3, -5, -10 



Acknowledgments 

Illuminating discussions with S. S. Ashwin are grate- 
fully acknowledged. 



[1] H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. 
Mod. Phys. 68, 1259 (1996). 

P. G. de Gennes, Rev. Mod. Phys. 71, S374 (1999). 
E. R. Nowak, J. B. Knight, E. Ben-Naim, H. M. Jaeger, 
and S. R. Nagel, Phys. Rev. E 57, 1971 (1998). 
S. F. Edwards and R. B. S. Oakeshott, Physica 157, 1080 
(1989). 

A. Mehta and S. F. Edwards, Physica A 157, 1091 
(1989). 

S. F. Edwards and C. C. Mounfield, Physica A 210, 290 
(1994). 

M. Schroter, S. Nagle, C. Radin, and H. L. Swinney, EPL, 
78, 44004 (2007). 

R. K. Bowles and I. Saika-Voivod, Phys. Rev. E 73, 
011503 (2006). 

S. S. Ashwin and R. K. Bowles, Phys. Rev. Lett. 102, 
235701 (2009). 



R. K. Bowles and S. S. Ashwin, arXiv:1102.0352 



There is no need for using two v-tile symbols. Both tiles 
have the same volume and energy. One cannot be ex- 
changed for the other in any given location. 
F. D. M. Haldane, Phys. Rev. Lett. 67, 937 (1991). 
Y.-S. Wu, Phys. Rev. Lett. 73, 922 (1994). 
S. B. Isakov, Phys. Rev. Lett. 73, 2150 (1994); Mod. 



Phys. Lett. B 8, 319 (1994). 

[15] D.-V. Anghel, J. Phys. A 40, F1013 (2007); Europhys. 
Lett. 87, 60009 (2009). 

[16] P. Lu, J. Vanasse, C. Piecuch, M. Karbach, and G. 
Muller, J. Phys. A 41, 265003 (2008). 

[17] D. Liu, P. Lu, G. Muller, and M. Karbach, Phys. Rev. E 
84, 021136 (2011). 

[18] P. Lu, D. Liu, G. Muller, and M. Karbach Condens. Mat- 
ter Phys. 15, 13001 (2012). 

[19] D. Liu, J. Vanasse, G. Muller, and M. Karbach, Phys. 
Rev. E 85, 011144 (2012). 

[20] The function S(V) can also be derived from |9j| via 
8 = -(dG/dT) p and V = (dG/dp) T , where G = 
— limjv->oo A r ~ 1 fcsTlnZ and the K m depend onp/T as 
noted. 

[21] Contact with quantities in Table I of [4] are readily estab- 
lished. If we set X = iff 1 , V = V/qa, S = S/k B , and 
w = wj, implying V = (2 + e 1 ^ = w 2 /(l + w), 

S = ln(l + w" 1 ) + V/X, and Z = (1 + w' 1 ), we obtain 
Y = V-XS = -llntl + itT 1 ) and, hence, Z = e"^. 

[22] Spatial correlations will be the focus of a separate study. 

[23] S. S. Ashwin, J. Blawzdziewicz, C. S. O'Hern, and M. D. 
Shattuck, Phys. Rev. E 85, 061307 (2012). 



